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Abstract 

This paper has a two-fold objective: determination of yearly momentum accumulation due to solar radiation Pressure 
and optimum reaction wheel sizing. The first objective is confronted while determining _propel ant 
the attitude control system over a spacecraft's lifetime. This, however, cannot be obtained from the daily momentum 
accumulation and treating that constant throughout the year, because the orientation of the solar arrays relative : to the 
spacecraft changes over a wide range in a year, particularly if the spacecraft has two arrays, one normal and the other 
off-normal to different extent at different times to the sun rays . The paper therefore first develops commands for the 
arrays for tracking the sun, the arrays articulated to earth-pointing spacecraft with two relational degrees of freedom 
and spacecraft in an arbitrary circular orbit. After developing expressions for solar radiauon torque due to one or both 
arrays arranged symmetrically or asymmetrically relative to the spacecraft bus, momentum accumulation over an 
orbit and then over a year are determined. The remainder of the paper is concerned with designing reaction wheel 
configurations. Four-, six-, and three-wheel configurations are considered, and for given torque and momentum 
requirements, their cant angles with the roll/yaw plane are optimized for minimum power consumption Finally 
their momentum and torque capacities are determined for one-wheel failure scenario, and six configuratio 
compared and contrasted. 


1. Introduction 

This paper is concerned with: a) determination of 
momentum accumulation due to solar radiation torque acting on 
an earth-observing spacecraft with sun-pointing solar arrays, and 
b) reaction wheel sizing and its pyramid configuration 
optimization for maximum momentum storage and minimum 
power consumption. These topics are classical; yet it seems 
there is no single reference in the published literature that treats 
this subject with sufficient comprehensiveness so that a control 
engineer, confronting this task, could accurately size the 
reaction wheels, select a wheel configuration, and estimate 
yearly propellant consumption for momentum dumping, all 
without extensive or expensive computer simulation. This 
paper, hopefully, fulfills that need. The contents of the paper 
and related previous contributions known to this author are 
summarized below. 

Solar arrays’ influence on spacecraft configuration, atutude 
control system, and mission operation is so far-reaching that a 
brief elaboration of this topic appears in order. For an earth- 
pointing three-axis stabilized spacecraft rotating once per orbit 
about the orbit normal, an attached solar array must at least 
have one relative rotational degree of freedom about the orbit 
normal so that the array can be held inertially Fixed and sun- 
pointing. Although economical, this arrangement becomes 
inadequate if the spacecraft's life span is more than several (say, 
six) months, because in this duration the earth moves around 
the sun in the ecliptic plane so much that the sun-rays deviate 
significantly away from the array normal and therefore a 
considerable power loss begins to occur. Consequently, for 
spacecraft with one year or longer life span, the solar array is 
accorded a second degree of freedom in the form of spacecraft 
yaw rotation. If the spacecraft at hand has only one solar array, 
this yaw rotation is of 180 degrees and may take place once in 
six months when the off-normality between the sun-rays and the 
array on one side of the orbit plane exceeds limits. The 180- 
degree rotation takes the solar array to the other side of the orbit 
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plane where the off-normality in the following six months will 
be smaller. Although to minimize cost and to gain simplicity, 
there are many spacecraft with one solar array (TOPEX for 
example), this configuration is asymmetric and might generate 
significant disturbance torques on the spacecraft arising from 
solar radiation pressure, gravity gradient, atmospheric drag, and 
thermal shocks at each sun-nse and sun-set; savings are 
therefore somewhat offset by stronger control torque 
requirements. To eliminate this asymmetry and/or to generate 
enough power for on-orbit needs, two solar arrays, one on each 
side of the orbit normal, are sometimes employed. Then, instead 
of once in six months, a yaw rotation from zero to 2k or from 
-n/2 to +k/2 takes place as continuously as one about the orbit 
normal. The corresponding sun-tracking commands were derived 
by McElvain (1961) 1 and Kalweit (1983) 2 . GPS satellites have 
opted for this approach. Although TOPEX satellite has one 
solar array, it also employs continuous yaw rotation instead of 
180° yaw rotation. Some missions are not interfered with by 
these persistent yaw rotations, but others are. For these latter 
situations, Kalweit 2 has determined best-fit minimum-power- 
loss, average yaw angles, constant over each half orbit. 
Nonetheless, persistent yaw rotation of a spacecraft is 
cumbersome because the torque and momentum capacity of the 
reaction wheels, usually employed for attitude control, must 
now accommodate the yaw rotation. A superior alternative 
appears to be, at least on the basis of technical merits if not 
cost, to bestow each array with two rotational degrees of 
freedom relative to the spacecraft, one about the orbit normal 
and the other about an axis in the orbit plane. Such is the 
spacecraft configuration considered in this paper; that is, a 
spacecraft with two solar arrays, arranged symmetrically on each 
side of the orbit normal and each array having two articulation 
degrees of freedom. Section 2 of the paper furnishes sun- 
tracking commands about the two just-mentioned axes. Explicit 
relationship is furnished between the so-called beta angle (also 
called flap angle) of the array and parameters such as earth's 
position in the ecliptic plane, the angle between the ecliptic and 
the equator planes, inclination of the spacecraft orbit, and its 
ascending node angle. 
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Turning our attention to solar radiation torque on a space 
vehicle, this arises from arrays as well as the vehicle bus. 
Moreover, of the two arrays, one may be normal to the sun 
while the other may be off-normal (thermal requirements may 
dictate so), and the bus and the arrays may cast shadow on each 
other at different times, changing the lit area, center of pressure, 
and moment arm from the vehicle mass center to the pressure 
center; see an example of shadowing in Evans (1964)-*. These 
complex effects are formulated and illustrated in Section 3. For 
typical spacecraft however, the torque contribution of the bus 
and the shadow effects are secondary; so by ignoring them, 
simple radiation torque expressions for arrays normal as well as 
off-normal to the sun are obtained which are used in Section 4 
to determine secular and cyclic momentum accumulation in the 
roll-yaw plane and about the pitch axis. Because of significant 
variation in the array's flap angle over a year, the corresponding 
momentum accumulation over each orbit changes considerably, 
particularly if the array is positioned off-normal to the sun. A 
simple expression of annual, secular momentum accumulation 
is therefore developed and illustrated in Section 4, and its 
dependence on the orbit inclination is investigated. From this 
result, yearly propellant consumption for momentum dumping 
is determined easily. DeBra and Cannon (196 1) 4 have also 
performed preliminary analysis along these lines. 

Section 5 of the paper is concerned with sizing reaction 
wheels and optimizing their pyramid configuration, keeping 
cost and redundancy in mind. Four-, six-, and three-wheel 
configurations with and without one wheel failure are 
considered. Optimum cant angles for these pyramid 
configurations for minimum power consumption and for given 
ratios between the roll, pitch, and yaw torque requirements are 
determined. Simple relationships are developed relating 
momentum or torque requirements about spacecraft axes to 
those about the wheel axes for all configurations with and 
without one wheel failure. These relationships then provide the 
required momentum and torque capacities of the wheels. The 
paper is finally concluded in Section 6. 

2. Commands for Sun-Tracking 
Coordinate Transformations 

In order to express sun-ray direction from the sun to the 
Earth in terms of solar arrays’ frames, and to develop pointing 
commands for the arrays for tracking the sun, the following 
angles, all anticlockwise positive unless stated otherwise, are 
introduced. The angle v, measured from the the first day of 
autumn (September 23), denotes the earth’s rotation around the 
sun in the ecliptic plane. The clockwise positive angle X ( = 
23.44 degrees) about the Vernal Equinox is the angle between 
the ecliptic and the equatorial plane. The angle Qn and i are, 
respectively, the usual ascending node angle of the spacecraft 
orbit in the equatorial plane and the orbit inclination angle from 
the equatorial plane. In this paper we will be concerned 
exclusively with circular spacecraft orbit. The local- vertical - 
local-horizontal (LVLH) frame £F C : 5Q Y c Zq at any point in 
the orbit locates the spacecraft mass center with X c along the 
velocity vector of the spacecraft, Z c along the local vertical 
from spacecraft to the earth, and Y c opposite to the orbit 
normal. To maintain the earth-pointing attitude, the spacecraft 
rotates clockwise about Y c -axis at the rate -coo (coo thus is a 
positive quantity and it equals the orbit rate of the spacecraft). 


The frame F c is the standard roll, pitch, yaw frame of a 
spacecraft with these three attitude angles zero. When the angles 
are nonzero, the spacecraft frame is denoted 3 r0 : Xo Yo Zo, as 
shown in Fig. 1. In this paper however, we assume that the 
spacecraft is controlled perfectly, and it always maintains its 
ideal LVLH orientation. 

We now define the orientation of the solar arrays relative to 
the frame or, equivalently, £F°. As stated in the Introduction 
and portrayed in Fig. 1, the two arrays turn relative to the 
spacecraft at the two-degree -of- freedom hinges Oj and O2. 
Considering +y-array first, its relative rotation is measured from 
the frame X10 Y\q Zjq which is parallel to the spacecraft frame 

In order to track the sun, the first rotation of the array is 
01y about the longitudinal axis y 10; this rotation annuls the 
clockwise rotation coot of the earth-pointing spacecraft measured 
from the ascending node line. The second rotation (often called 
beta angle), denoted here 0 1 z , takes place about the once- 
displaced Zio-axis. We thus arrive at the array-fixed frame fF 1 : 
Xi Y\ Zi, with the array in the YjZi plane and its outward 
normal along Xi. Note that when the array is normal to the 
sun, the sun vector £. from the sun to the earth is opposite to 
the array normal X\. The transformation between the frame £F® 
and the array frame ^F 1 for the sequence 0j y , 0i z is 


' Xo 


C01y C01z 

— C0 ly S0 1 z 

S01y 


Xi 

Y 0 

= 

S0 lz 

C0tz 

0 


Y, 

L Zo J 


-S01yC01z 

S0ly S0lz 

c 0 1 y j 


L A J 


where c(«) = cos(-) and s(») = sin(«). The initial orientation of the 
-y-array is the same as that of the +y-array, and the rotations of 
the -y array are conveniently measured relative to the frame 
X20Y20Z20 with Y20 opposite to Yo and Z20 opposite to Zo. 
The frame X20Y20Z20 * s selected such that the solar cell face 
of both arrays are on the same side. The first rotation 02 y about 
the Y20-axis (Fig- 1) nullifies the orbit rotation, and the 
rotation 02 z about the edge Z2 brings the array to the desired 
normal orientation relative to the sun. The transformation 
matrix from the spacecraft frame £F° to the array-fixed frame 
& 2 : X 2 Y 2 Z 2 is 


Xo ^ 


C02zC02y -S02zC02y S02y 


' x 2 ’ 

Yo 

= 

_s ®2z “C 02z 0 


Y 2 

Zo _ 


C02*S02y -S0 2z S0 2 y -C0 2y 


. A - 


Commands For Sun-Tracking 

Because the solar arrays are hinged to the spacecraft, it is helpful 
to express the unit vector £ from the Sun to the Earth in the 
spacecraft frame £F°, or equivalently, in the orbit frame £F C . Let 
S c i, S c 2> and S c 3 be the components of 3. in £F C . Then, using 
coordinate transformations involving the angles v, X, On, i, 
and coot, defined above, these components are found to be 
S c i = coot (ci(-s£2NCV + cQfsfcXsv)+sisXsv} 
-soot(cftNCV + cXsv) 

S C 2 = si(-sftNCV + cftNctevJ-cisXsv 

S c 3 = -scoot {ci(-sI2Ncv + cQivjcA.sv)+sisXsv) 

-coqUcQncv + sQn cXsv) (3) 
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Fig. LA Spacecraft with +Y and -Y Arrays, their Frames, and 
Articulation Degrees of Freedom 

+ Y-Array Commands 

Because the first rotation 0i y about the longitudinal axis 
y 10 annuls the once-per-orbit rotation of the earth-pointing 
spacecraft, it may be clear that, ideally, 

0 1 y (0 = (OQt + 010 W 


Eq. (8) indicates that, as the earth rotates around the sun (0 < v 
< 2k), the angle 0i z varies sinusoidally with a certain 
amplitude; this is illustrated below in three examples. 

Example l : = 0, v = 7112 

The correctness of Eq. (7) and Eq. (8) can be illustrated by 
considering v = k/2 and Qn = 0° for which 

010 = tan -1 0/[ — c(X — i)] = 7t (9a) 

0 1 z = ^ — i (9b) 

If we further assume that the satellite orbit lies in the ecliptic 
plane, which means the orbit inclination i equals X, we will 
have 0i2 = 0. In Eq. (9a), the choice 0 to = k and not zero is 
selected to ensure that the solar cell face of the array, not its 
back side, is towards the sun. Moreover, Eq. (9a) and 0i z = 0 
together imply that in order to be normal to the sun rays at t = 
0, the array must be in the plane Y c Z c . Physically, the 
conclusion 0 1O = ±re and 0 lz = 0 is seen to be valid in Fig. 2 for 
the parameters v = k/2 , = 0, and i = X at t = 0. 

Example 2 ; - 0, and arbitrary v 

For ftjsj = 0, Eq. (7) and Eq. (8) yield 


where 0io is 0i y (t=O). Note that 0 i y (t) is anticlockwise 
positive, whereas the negative sign of coot has already been 
accounted for in (3), so coot in Eq. (4) is positive: coo^O. To 
determine 0io and the second rotation 0] Z , we observe that 
when the array is normal to the sun, the incoming sun vector S. 
is opposite to the outgoing array normal X i : 

= [X] Y, Zi ] T = [-1 0 0] T (5) 

where the superscript T means the transpose of the column 
vector. Substituting Eq. (5) in Eq. (1) the sun-ray vector £ in 
the LVLH frame is found to be 


' S cl l 

cfO 

-c0 ly c6] Z 

S c 2 

= 

-S01z 

. S C 3 . 


S0 ly C0 l z 


010 = tan -1 


-cv 


-c (X-i) sv ’ 


9lz = sin 1 


[ sinv sin (X-i) ] 

(lOa.b) 


Eq. (10b) states that, inasmuch as (X-i) is fixed, the angle 0i z 
will change periodically as the earth moves around the sun in 
one year, 0 < v < 27t; the extremes of 0i z will be (X-i) when v 
= k/2, and (i-X) when v = 2>k/2. This is illustrated in Fig. 3 for 
i = 28.5° and 43.5°. Now considering the satellite orbit in the 
ecliptic plane (X=i), Eqs. (10) simplify and, consistent with 
Example 1, furnish 


0 ! 0 = - k /2 - v (11a) 

0 1 y (t) — co 0 t — (k/2 + v) (lib) 

For v=0, +y array's rotation 0io = -k/2 about the axis Y c , 
measured from its datum orientation in the Y C Z C plane, is 
illustrated in Fig. 2. Also, for Qjq = 0 and i = 28.5° the linear 
relationship (11a) is seen to be true in Fig. 3. 


The unknown initial angle 0io is determined by substituting 
Eq. (4) in Eq. (6), yielding: 

_ 1 -[ cQn cv + sQn cX sv] 

10 - tan c - cv + cftjsf cX sv) + si sX sv ] 

(7) 

where the negative signs in the numerator and denominator are 
retained so as to arrive at a unique value of 0 1O within the range 
-K < 0] O - k; this will ensure that the solar cell face of the 
array looks at the sun. The array’s inclination angle 0 1Z is 
obtained from the second components of Eqs. (3) and (6): 

012 = sin -1 (si (sX2n cv * C ^N sv ) + ci sX sv] (8) 


Z VE Ae 



WE 


Fig. 2. +y Solar Array Orientations at t=0 for Examples l and 2 
(VE = Vernal Equinox) 
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(A) j - 28 5°, \ m 23 44“ (8) j - 43 5 a , X - 23 44° 



f7g. 5. Initial Angle 9jo about the pitch axis yo and inclination angle &i z about the short edge Zj-cvas of +y array versus angle v 


Table L Extrema of the tilt angle Gj 2 


Earth’s Position 
in the 

Ecliptic Plane (v.deg) 

Ascending Node Angle 
of the Spacecraft's Orbit 
(Qn. deg) 

Extreme 
Value 
of 8iz 

90° 

0 0 1 

X-i 

90° 

180° 

X+i 

-90° (or 270°) 

0 0 

i-X 

-90° (or 270°) 

180 0 

-X-i 


Example 3: Arbitrary and v 

Fig. 3 illustrates the variation of 0io and 0iz as a function 
of the earth’s position in the ecliptic plane ( 0 < v < 2n ) for i 
= 28.5° and 43.5° and eleven values of the ascending node angle 
Qfsj. The apparent discontinuities in 9 1 o curves at ±rc are 
inconsequential because the angle 0io has a range of 2 n and +rc 
= -k. Furthermore, from Eq. (8), we infer that 


It is instructive to compare the commands developed above 
with those developed by McElvain 1 and Kalweit^. 

3, Solar Radiation Torque 
Radiation Torque on a Solar Array 

For momentum accumulation study, the radiation torque g 

at the spacecraft mass center is required. Referring to Fig. 1, let 
t j (j = 1,2) be the vector from the reference origin O to the solar 

array hinge Oj (j = 1,2), and P p j (j = 1,2) the vector from the 

hinge Oj to the pressure (or geometric) center of the array. 
Additionally, let rc be the vector from the reference origin to the 
vehicle mass center. Denote the moment arm vector of the solar 
radiation force on the j -array as £ pj . Then, 

£pj = -Ic + bj + Coj ^pj (j=U) ( 16 ) 


for Q N = 7t: 9 1Z = sin 1 [ s (X+i) sv ] (12) 


which is a counterpart of Eq. (10b) in Example 2. Clearly, the 
extremes of 0i z are 


f X+i @ v = n/2 
\-(X+i) @ v = 3rc/2 


(13) 


These, as well as the sinusoidal variation of the angle 0iz 
versus v, are illustrated in Fig. 3. Also, see Table 1. 

-Y Array Commands 

Regarding the -y -array, because the angles 02y and 02 z are 
defined about the axes Y 20 and Z 20 which are respectively 
opposite to the +y-array axes Yjo and Z 10 , it is clear that for 
keeping the array normal to the sun-rays 

02y = — 0 1 y (14a) 

02z = -8iz (14b) 


Recalling Eq. (4), the desired 02y(t) will therefore be 

02y(t) = -coot -010 (15) 


where Q>j is the transformation matrix defined by Eq. (1) and 
Eq. (2) for j=1.2. Following Reference 5, Section 8.3, the 
radiation torque experienced by a spacecraft about its mass center 
owing to the sun rays off-normal to the jth-array, is 


Jj = pAj cctj £pj X [ (o t + c rd ) £ + 2 (Orf/3 + o n cap flAj ] (j=l,2) 

(17) 

where p denotes the radiation pressure on a totally absorbing 
normal surface, Aj is the array s area, cctj = cos ctj and ctj is the 
angle between the inward normal unit vector and the sun 
vector £, a a is the absorptivity coefficient of the surface under 
consideration, c r d is the diffused reflectivity coefficient, and a rs 
the specular reflectivity coefficient. The notation x in Eq. (17) 
transforms a column vector to a 3x3 skew-symmetric matrix, as 
defined by Eq. (13), Section B.2, Reference 5. 

Array Normal to the Sun 

When the sun rays are normal to the array, caj = 0, and 
IlAj = £ and Eq. (17) simplifies to 
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where 


gj = pAj CT A £pj X £ 

(18) 

<7 A = 1 + <Trs 2a rc j/3 

(19) 


The three components of the radiation torque, Eq. (18), on the 
spacecraft in an arbitrary orbit, with the array normal to the sun, 
will now be shown explicitly. For concreteness, we assume that 
the spacecraft mission is such that the power and thermal 
requirements allow the +y -array to be normal to the sun. To 
express the vector £ p i component- wise, note that the length of 
the A-frame along the pitch axis is ‘a’ (Fig. 1), and the rotation 
01 Z takes place about the Zi-axis of the array without involving 
the A-frame. Because in the X} YjZi frame, the pressure center 
of the array is at a distance Z/2 along the yi-axis from the 
transverse edge of the array (Fig. 1), the vector £oi P p i in Eq. 

(16) can be calculated easily using Eq. (1). Furthermore, we 
assume that when the arrays are in the yozo plane (0 lz = 0 = 
02 Z ), the vectors ^ and b] are in that plane (Fig. 1). Therefore, 
in the LVLH frame cF c : 


- Ic + hi — [ 0 “fcy + hly TcZ+^lZ (20) 


The vector £ p i, Eq. (16), will therefore be 


£pl = 


-b S0 1 Z C0 1 y 
-Tcy + bi y + a +b C0 1 z 
-r C z + bi z +b s0iz s0iy 


b* 1/2 

( 21 ) 


For analytical convenience, define 

biy® = -icy + bly + 3. = — r cz + bi z (22) 

Using (21), (22)in Eq. (18) the following components of the 
torque g ^ in the LVLH frame are obtained: 


gj = P*! a A 


b® S0 lz +(b+b® C0 lz ) S0]y 
“ bj z C 0 lz C0 1 y 
(b+bf y c0 lz )c0ty 


(23) 


While the angle 0i y (t) varies linearly at the rate co 0 , changing 
by 2 n over one orbital period, the angle 0i z is virtually 
constant in that period. Eq. (23) therefore indicates that, in the 
LVLH frame, the roll (x-) torque comprises a constant and a 
cyclic torque, the pitch and yaw torques are cyclic, and the cyclic 
yaw torque is in quadrature with the cyclic component of the 
roll torque. (Also, see Ref. 4.) 

-Y- Array Off-Normal To The Sun 

When the spacecraft at hand has two arrays, power and 
thermal requirements might dictate one array to be normal to the 


sun rays and the other array off-normal. Fig. 4 depicts one such 
posture for the spacecraft in Fig. 1, with +Y-array normal and - 
Y -array off-normal. 

To keep the mass center of the -y-array on the rotational 
axis yoi the A-frame is turned about the Z2-axis by an angle 
^2z,a ( a negative 02z,a is shown in Fig. 4); also, the array’s 
normal -to- the- sun orientation angle is denoted 02 z » whereas the 
off-normal orientation angle is denoted 02 z '- Under these 
circumstances, the transformation matrix £o2> defined by Eq. 
(2), is altered to 0)2, A for the A-frame replacing 02 z with 
02z,a» and to £o 2 for the -y-array substituting 02 z ' for 02 z 

Thus the pressure center vector £02 Pp2 i n Eq. (16) is revised 
to 



~S02z,a c 02y 


~S02z c 02y 

Qo 2 P P 2 = = a 

“0^2z,a 

+ b 

~O02z 

--s02z,aS0 2 y - 


S02z s ®2y - 


(24) 


To evaluate the torque gj 0=2), we note that the vector - r ^ + 
£t2 for the -y-array in the frame £F C is, following Eq. (20) 

-Ic + 122 = [0 -r cy + b2y -r cz + b2 z ] T (25) 

where, to symmetrize the mass distribution 

b2 y = -bi y b2 z = bi z (26) 

Using the equation 02y = -©ly and Eqs. (24), (25), and (26), the 
vector £ p j (j = 2), Eq. (16), is found to be 



-b s 00 ] y 

gre 

c p2x 


— r C y -bl y — b c 

m 

Cp2y 


~ *"cz ^“^lz S01 y 


_ 0p2z 


%>2 = 

where the lengths b c and b s are defined as 


be — a C02 z ,a + b C02 z 
b s = a S02z,a + b S02z' 


gx 


(27) 


(28) 


The vector £ and the inward normal qa for the -y-anay in 
the off-normal frame (F 2A : X2A ^2 A z 2A (Fig- 4) are 


= [— c§2z -sS* o ] T , Oa^ 2 * = [ -1 0 0 ] T (29) 


where 822 is the off-normal angle, defined by 

62* = 0^ -0^' (30) 

and equal to the angle ct2 between the vectors £ and ha for the 
array. Therefore the vector [ • ] in (17) the spacecraft frame is 
found to be 
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Fig . 4. A Spacecraft With +Y Array Normal to the Sun and -Y 

Array Off -Normal 


_ 

[ (<T a + °id) S. + 2 (O r d/3 + O rs C(X2) Ha ] ° = 

[ -O cx C02y Ocy -°cx S&2y 1 T (31) 

where 

Ocx = 0x2 C02z ' - Oy2 S02z' 

O cy = O x2 S02z ' + Oy2 C02z' 

o x2 = (1 + o ra ) C& 2 , + 2o rd /3 

o y2 = (l -CTjsS^ (32) 


Using the vector £p2 ' n Eq. (27) and the vector Eq. (31) the 
torque §2 is found to be 

g 2 3rC = P A 2 cS2z 


- (— r cz + b lz ) o cy - (b s o cy + (r cy + b ly + bj o c * } s0 ^ y 

— (— 1 Tcz + biz) Ocx ) C0 X y 
- (bs Ocy + (r C y + bl y + be) Ocx )c0 1 y 


(33) 


which exhibits the same attributes as those exhibited by _g ^ in 
Eq. (23). 

Example 1 : Off-Normal Angle &2i - 0 and A-Frame Angle $ 2 i a = 0 
In this circumstance, the +y- and -y-array become parallel 

because now &2z' = -0iz» and the torque g2^ C i Eq. (33), 
simplifies to 

\ b®se u + (-b + b® ce u )se ly 1 


= pA 2 Oa 


-b 2z c9uc0i y 
(-b + b® c0 u )c0 ly 


(34) 


Example 2: Both Arrays Normal to the Sun: Resultani Torque 

Adding Eq. (34) to Eq. (23), and assuming that the arrays 
are identical in geometry as well as in optical surface properties 

so that A s = A] = A2, the resultant torque g equal to E 1 + 
g 2 ) is found to be 


g& c = 2pA s o A 


(-r C z + bu)s0iz -r cy C0 ljt s0 ly 

( r cz — ^1 z) tz C® 1 y 
— Tcy C0 l z C0 1 y 


(36) 


Mass asymmetries in the spacecraft generate the components r cy 
and r cz ; how each component contributes to the radiation torque 
on the spacecraft is seen clearly in Eq. (36). We also observe 
that, because the two arrays are parallel and geometrically 
identical, the A-frame length a and the array’s half-length b do 
not contribute to the total torque, only the hinge location 
referenced from the vehicle mass center matters. 

Example 3: +y- Array Normal and -y - Array Off-Normal: Resultant 
Torque 

The resultant torque g :f c is now obtained by adding Eq. 

(23) with Eq. (33). To build a simple expression for jL^ c , 
observe that, in the £F C frame, ignoring the variation in the 
angle 0| Z over one orbit period, the constant part g c of the x- 
torque (roll torque) is found to be 

g c = pA s b) z 0 [a A s9 1 z - a cy c52 z ) (37a) 

and the amplitude g a of the sinusoidally varying part g a s9i y of 
the x- torque is 

ga = pA s [<3 a (b + biy® C0| z ) - c82z (b s G C y + 

(r cy + b| y + b c ) a cx ) ] (37b) 

The z-component (yaw torque) has no constant part, and its 
sinusoidal variation is in quadrature with the roll component. 
Lastly, the y-component (pitch torque) also varies sinusoidally, 
with its amplitude gb equal to 

gb = -pA s biz® [°A c ®lz + °cx c &2z] (37c) 

Thus, the total radiation torque g - 7C acting on the spacecraft at 
its mass center, in the spacecraft frame, is 


gc + ga S01y 


g 


^C 


gb C01y 
ga C 01 y 


(37d) 


where 

t>2y 0 = -r cy _ 6] y -a b 2 z ® = b] z e (35) 


The constant roll torque g c in the rotating frame & c becomes a 
periodic torque in the orbit plane in a non-rotating frame; 
conversely, the periodic roll-yaw torque in 3 :c becomes a 
constant torque in the orbit plane in an inertial frame. 
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Meanwhile, the periodic pitch torque remains periodic even in a 
non -rotating frame because the rotation (Dot takes place about 
the pitch axis. 

Radiation Torque on Spacecraft Bus 

Usually the surface area of spacecraft bus is much smaller 
than that of the deployed solar arrays, and therefore the radiation 
torque caused by the bus is an order of magnitude smaller than 
that caused by the arrays. For this reason and to conserve space 
the analytical details of derivation of the radiation torque on a 
bus will not be given here. 

Torque On An Array Shadowed By the Spacecraft Bus 
Fig. 5 portrays a spacecraft with a cylindrical bus and a 
solar array, the bus casting a shadow on the array. The lit semi- 
cylindrical surface is ABC. The shadow boundary on the array 
may utmost consist of three segments: an arc cast by the lit 
end-face of the cylinder and two straight lines, enclosing the arc, 
cast by the boundary rulings at A and C. The shadow in Fig. 5, 
though, has two segments only because the shadow of the 
ruling at A falls off the array. The shadow boundary is 
determined as follows 

It may be clear that at any time only one array at the 
most will be shadowed. Knowing the angle (5 between the sun- 
ray unit vector S. and nadir unit vector £ 3 , we first determine 
which end-face is lit: 

if fj < rc/ 2 , the earth-pointing end-face is dark, and the 
opposite end-face is lit; 

if p > 7 t/ 2 ,the earth -pointing end-face is lit, and the 
opposite end-face is dark 

and the lit semi-cylindrical surface is then identified. Next, let r 
be the vector from the reference origin O to any point P on the 
edge of the fully-lit face of the bus or along the two boundary 
rulings of the lit semi-cylinder. The vector from the hinge Oj of 
jth array (j = 1,2) to the point P is then -bj + r. The two 

components of this vector in the array plane yjZj, which is not 

necessarily normal to the sun, are 

(-bj + r) . ij and (-bj + l) • kj (38) 

where jj and kj are the unit vectors along yj- and Zj-axis of the 
array, respectively. If these components fall on the array, that 
is, 

a < (-bj + D ■ ij (39a) 

- w /2 < (-bj + r) . kj < w /2 0 = 1*2) (39b) 

where w is the width of the array along the Zj axis, then the 
array j is clearly shadowed. In that event, as the vector r moves 
along the edge of the lit end-face and the boundary rulings, the 
above two components will delineate the shadow boundary. 

The radiation torque due to the array is still given by Eq. 
(17). The area Aj represents the lit area of the array, and the 

vector Ppj in the moment arm £pj, Eq. (16), is the vector from 
Oj to the instantaneous center of pressure of the shadowed array. 

The instantaneous lit area and the vector P pj are calculated 
numerically. 


•-QWAAOS 
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Fig. 5 Cylindrical Spacecraft Bus Casting Shadow on the Solar 

Array 


ALL TORQUES IN FT-LB 


SOLARTORQUE SHAOOW 


-000025 PERIOD 




- = WITH SHADOW I 

- NO SHADOW I 

Fig. 6. Comparison of the total radiation torque on the 
spacecraft with and without considering the shadow; Orbital 
parameters : v - 90°, Qpj = ISO 0 , i = 20°, corresponding solar 
arrays' inclinations: 0j z = -&2z = 43.5° (arrays parallel and both 
normal to the sun) 


9 


HH/VT-032 




3001 



Yj-AXIS (inches) 


Fig . 7. Shadowing of -y array: Orbital Parameters and solar 
arrays' inclinations same as those in Fig . 6 

Illustration 

For parameters not recorded here because of space 
limitations. Fig. 6 furnishes total radiation torque (two arrays 
plus cylinder) with and without considering the shadow of the 
bus on the array, for ^=0 and both arrays normal to the sun 
rays. When the shadow is ignored, we observe that, in 
conformance with Eq. (36), the yaw torque is zero and the roll 
torque a constant over a few orbits. Slight oscillations in the 
roll torque are present because of the cylindrical bus. The 
variation of the pitch torque in Fig. 6 conforms with the y- 
component of the solar radiation torque in Eq. (36). Regarding 
the effect of shadow on the array, we observe that while the 
shadow does not alter y-torque, the trough of the x-torque 
plummets from -0.29E-3 ft.lb to -0.38E-3 ft.lb during the 
shadow period. A small (0.1 E-4 ft.lb amplitude) cyclic z-torque 
also arises during the shadow period. For the parameters under 
consideration, it turns out that the -y-array is shadowed, +y- 
array is noL Traversal of the cylinder s shadow on the -y-array 
is shown in Fig. 7. The growth of the shadow area on the -y- 
array is shown in Fig. 8 for three parametric sets, including the 
set for Fig. 6. Fig. 8a shows that, for example, when 02 z = - 
43.5°, the maximum shadow area equals 88 sq.ft [(shadow area) 
/ (array area) = 88/202.85 = 0.43] and the array remains partially 
shadowed for nearly 6 hours (one-fourth of the orbital period). 
We nonetheless also observe that the shadow occupies 43% area 
of one array only briefly (several minutes), occupying 
progressively smaller area before and after the maximum shadow 
epoch. As the shadow traverses, the pressure center of the lit 
portion travels also; the loci of the instantaneous pressure 
center, from the moment the shadow enters the array till the 
moment it leaves, are shown in Fig. 8b for three sets of 
parameters. As expected, these loci are closed curves, and the 
bigger the angle I02 z h the wider the loci. The shadows occur 
around different orbit angles coot for different sets of orbital 
parameters v and On (even if the orbit inclination i is the same 
for these orbits). For plotting convenience, however, the 
maximum shadow epoch is shown to be the same (t = 18 hrs) 
in Fig. 8a for all three sets of parameters. In the preceding 
results, the center of mass vector ^ is zero. However, as the 
arrays orientations change relative to the spacecraft bus while 




Fig. 8 . (A) Shadow area on the -y-array versus time; (B) Locus 
of the pressure center of the -y-array as shadow travels on the 
array; orbital parameters for 0 i z = -43.5 ° the same as those in 
FigJ; for $ 2 Z = -30°, the parameters are v = 30° and = 
140°; for = -22°, v = 80°, = 80°, and i = 20° 

tracking the sun, £c may vary depending on the arrays' 
arrangement. In the case at hand, this variation is found to be 
small, and it was not ignored in the computations. 

The preceding results show that the spacecraft bus shadow on 
the array does not change radiation torque significantly. Also, 
based on the results not shown here due to space limitations and 
as stated before, the radiation torque on the bus is found to be an 
order of magnitude less than, and so negligible compared to, 
that on the arrays 

4. Momentum Accumulation Due to Solar 
Radiation Torque 

In the following analysis we will determine the momentum 
accumulation due to the arrays’ torque only. It is instructive to 
compare the following development with that by McElvain 1 . 

Short Term Accumulation 

Let lib be the inertial angular momentum vector of the 
earth-pointing spacecraft under consideration. If H x , H y , H z are 
the components of Kb in the body-fixed frame, then, in the 
presence of the radiation torque, (37d), they will be governed by 
the following three equations, expressed in the S 70 : 
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H x - ©0 H z 

= gc+ gaS01y 

(40a) 

k, 

= gb c ^ly 

(40b) 

H z + ©0 H x 

= ga c ^ly 

(40c) 


As is well known, the roll (H x ) and yaw (H^ momentums are 
gyrically coupled, whereas the pitch momentum (H y ) is 
independent of roll and yaw. Examining the definitions of g a , 
gb> gc in Eqs. (37 a, b, c), we observe that they are functions 
of, among other things, the inclination angle 0j z and the off- 
normal angle 52 z * Because these angles depend on the earth’s 
motion around the sun (the angle v), they are held constant in 
flight over a few days, and therefore, within this period, the 
quantities g a , gb> gc can be regarded as constant. Moreover, 
assuming that time t can be measured such that 0i y (t=0) = 0, 
instead of 0io according to Eq. (7), 0i y (t) in Eq. (40) can be 
replaced with coot. With these two assumptions, Eq. (40a) and 
(40c) can be integrated together, yielding (i 2 = -1) 

H x + iH z = -i (gc/coo) (l-e - ' 0 ^ 1 ) + i t g a e- ,c °0t (41) 

The integration of Eq. (40b), on the other hand, yields the pitch 
momentum H y as 

H y (t) = (gb/©o) sin coot + H y o (42) 

Because a constant torque gc in the orbiting frame is a cyclic 
torque in the inertial frame, whereas a cyclic torque g a in the 
orbit frame is a constant torque in the inertial frame, it is 
natural to find in Eq. (41) a cyclic variation in (H x + iH z ) 
owing to gc and a secular growth owing to g a . At any arbitrary 
instant, the magnitude of the secular term in H x +iH z is tg a , 
showing a linear growth in the spacecraft momentum in the 
roll-yaw plane. The pitch momentum H y (t), Eq. (42), varies 
cyclically, with the amplitude equal to gb/wo. Th e constant 
Hyo in Eq. (42) equals H y (0) and it embodies the spacecraft’s y- 
momentum arising from its once-per-orbit equilibrium rotation. 

Yearly Accumulation 

It was just shown that the secular momentum in the roll- 
yaw plane at the end of one orbit (wot = 2 n) is (H x + ify = i 
2 tc ga/wo* where the coefficient g a , Eq. (37b), is a function of 
the solar array’s inclination angle 0i z for an array normal to the 
sun, and of the inclination angles 0'2 Z and off-normal angle $2z 
for an off-normal array. Fig. 3 shows that the angle 0i z (or 02 Z 
which is equal to - 01 *) varies as a function of the angle v (the 
earth’s motion around the sun); therefore in order to calculate 
the yearly momentum accumulation, this variation must be 
considered. To separate the linearly varying angle v(t), 0<v<2rc, 
from other constant orbital elements X, Qn> and i in the 
definition of 0i z , Eq. (8d) is rewritten thus: 

0 lz = sin“ l [ (si s^n) cv + (" si C ^N cX + ci s ^) = 

sin -1 [Ao sin (vo + vt)] (43) 

with the amplitude Ao and the phase angle vq defined by: 


A 0 = (A f 2 + V) 1/2 ^ 1; v 0 = tan 1 [AJA C ] (44) 

A s = si sQn A^ = -si c^n cX + ci sX (45) 

Strictly speaking, because of orbit regression the ascending node 
angle Qjq is not constant, but this variation is ignored here. 

Spacecraft With + Y -Array Only; 

Array Normal to the Sun 

When the +y-array is normal to the sun, the coefficient g a 
given by Eq. (37b) simplifies to 

g« = pA s OA(b + bi y ® c0iz) (46) 

ignoring the second term in [•] in Eq. (37b) because it pertains 
to the -y -array. The magnitude of the secular momentum 
accumulated at any instant t then becomes 

IH* + iH z l = pA s a A (b + bi y ® c0j z ) t (47) 

The yearly momentum is obtained by integrating (47) over the 
annual variation of the angle 0i z . Let n y be the current orbit 
under consideration: 

n y = t/xo , to = 2n/wo (48) 

where to is spacecraft orbit period, and let N y be the total 
number of spacecraft orbits in one year. Then the yearly 
accumulation will be 

Ny 

Hyr = (pA s a A to) J (b + biy® c0iz) dn y (49) 

0 

Because -nfl < 0i z < rc/2, c0i z will always be positive. With 
the aid of Eq. (43) and recalling that b is constant (= Lfl), and 
treating bi y ®, Eq. (22a), as constant, Eq. (49) transforms to 

Ny 

Hyr = (pA s a A To) [bNy + biy® J ( 1-Ao 2 sin 2 

0 

(vo + vxo n y ) } 1/2 dn y ] (50) 

(The assumption bi y a constant may not be always valid, 

because the quantity r cy involved in the definition of biy may 
change due to 0i z , but this change can be minimized by 
keeping the array’s mass center on the axis of rotation.) Now, 
while n y changes from 0 to N y , the argument of sin 2 in Eq, 
(50) varies from vo to vo+2rc; therefore, the integral in (50) is 
an elliptic integral of the second kind. Measuring the time t 
such that the phase angle vo = 0 and capitalizing on the 
symmetry properties of the elliptic integral at hand, after 
integration is found to be 

Hyr = (pA s <J A To) [bNy + 4 (bi y ®Mo) E (tc/ 2 , Ao)] (51) 

where E(tc/2, Ao) is a complete elliptic integral given by the 
series: 


11 


HH/VT-032 



H z (ft.lb.s) 
30 

* 6tl 




(52) 


Betas uc the number of the spacecraft orbits in one year is 


Ny = 2 h/(vt 0 ) (53) 

the yearly momentum accumulation Hy T simplifies to 

Hyr = Ny TO pA s a A [b + biy® (1 - Aq 2/ 2 2 - ...) ] (54) 

In the series (52) or (54), as many terms are retained as are 
necessary to evaluate the sum up to a desired accuracy. 

Spacecraft With +Y and ~Y Arrays . 

a) Both Arrays Normal to the Sun 

The coefficient g a , Eq. (37b), now simplifies to 


4rc 


23ji/4 


1 5ti/4 


24 
18 

12 ; 2K 
7kJ4 6 


1 7tu'4 


9n/4 


kJ 4 


i — - — 0-^ ^ iv H x (ft.lb.s) 

-30 _24 -18 - 12 , -6 0 6 ^ 12 , J 8 _ 2 * 30 


1 1 it/2 ‘ 7 tl/2 ' 3ji/2 


13tc/4 


-72- n 


-IB- 


21 ti/4 


-24 
-30 : 


Jt/2 571/2 9 ti/2 

' 3ti/4 

' 11 71/4 

19ti/4 


5k 


g a - -2pA s C?a r cy c ^lz (55) 

This is corroborated by the coefficient of s0i y in the first 
element of the vector equation (36). Comparing (55) with (46), 
the yearly momentum accumulation in the present case can be 
written down following Eq. (54): 

Hyr = N y to pA s <J A [ -2r cy (1- A<) 2 /2 2 - ...)] (56) 

h ) +Y-Array Normal , -Y-Array Off-Normal to the Sun 

In this case, the coefficient g a is given by the full-length 
equation (37b). The definition of the coefficients b c , b s [Eq. 
(28)], and a cx , a cy [Eq. (32] reveal the presence of the products 
of the trigonometric functions of the angles 0'2z and 02z,a* The 
analysis therefore seems intractable, and developing a closed- 
form expression for infeasible; consequently, numerical 
integration of the equation IH X + iH z l = t g a over one year 
seems inevitable. 

Yearly Propellant Consumption 

Annual propellant consumption, W p , should not be based 
on the momentum accumulation over one orbit and treating 
that constant for entire year, for the daily momentum varies 
significantly over one year (illustrated later in Fig. 10). Instead, 
knowing H yr from Eq. (54), Eq. (56), or otherwise, and 
knowing the specific impulse I sp of the propellant under 
consideration and moment arm lj of the thrusters from the 
vehicle mass center, W p is obtained from 

Wp = H yr /(I sp *j) (57) 

Illustrations 

Example L Spacecraft With Two Arrays: Momentum Accumulation 
Over Three Orbits 

Fig. 9 illustrates momentum accumulation in the roll-yaw 
plane, Eq. (41), with the orbit angle (Oot as parameter. In 


Fig. 9. Momentum accumulation in roll-yaw plane, due to the 
two arrays only: + Y array normal to the sun, 0} 2 = -45 -Y 
array off-normal, &2z' - 25°, &2z = 20°, -Y array A-frame angle 
&2z,a = -57° 

particular, it applies to the spacecraft configuration shown in 
Fig. 4, with +y-array normal to the sun and -y-array off-normal 
at an angle Fig. 9 shows the linear, radial growth of (H x + 

iH 2 ) — the term itg a e _,0) O l in Eq. (41) — superimposed upon a 
cyclic variation with an amplitude of gJu>o. 

Example 2. Momentum Accumulation As a Function of Earth’s 
Position in the Ecliptic Plane 

In order to obtain complete dependence of the daily 
momentum accumulation over a one year period, we next 
illustrate in Fig. 10 the roll-yaw momentum at coot = 2n and 
the pitch momentum amplitude gb/coo as a function of the 
angle v: 0 < v < 2 n y with the off-normal angle 5 2z as a 
parameter. Figs. 10a and 10b show that, for off-normal angle 
52 z = 0, the roll-yaw and pitch momentum amplitude varies 
periodically with v, with half-year period; the corresponding 
ratio H y / IH X + iH z l remains constant for entire year. When the 
off-normal angle §2z of the -y-array is introduced, the half- 
yearly periodic variation of the roll-yaw momentum disappears 
and, instead, it begins to vary asymmetrically with v. Also, 
the peaks and valleys of the roll-yaw momentum grow further 
apart as §2z increases. The half-year periodic variation of the 
pitch momentum amplitude alters little with the angle 52 z . Fig. 
10c illustrates the ratio H y / IH X + LH Z I; the variations in this 
ratio with v become more pronounced as §2z increases. This 
ratio is of interest because it helps decide the cant angle of the 
reaction wheel configurations considered in the next section. 

Example 3. Yearly Momentum Accumulation 

When the off-normal angle 6^ equals zero (that is, when 
both arrays are normal to the sun) the yearly momentum can be 
obtained analytically using Eq. (56), and when 8 2z * 0, the 
yearly momentum is obtained by numerical integration of the 
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area under the curves IH X + iH r l in Fig. 10a for the entire range 
0 < v < 2k. Fig. 1 1 depicts the yearly momentum against 5 2z 
for the range 0 < 8 ^ < 20 ° for different orbit inclination angles, 
keeping Q N at 170°. For 8 ^ = 0, the analytic prediction was 
compared and found identical with the numerical results. The 
number of terms that must be retained in the infinite series in 
Eq. (56) increases with the orbit inclination i, and yearly 
momentum diminishes a little. The yearly momentum, 
however, increases significantly with the off-normal angle 822 , 
as seen in Fig. 1 1 . 


O a * 0.75, CJrs-0.25, Ord = 0 0 ' A s1 = 202.85 ft 2 
i = 28. 5 a , On = 170.0, X * 23.5° 



15 90 135 160 225 270 315 360 

v (DEG) 


Fig . 10. Variation of roll-yaw and pitch momentum with the 
earths position (the angle v) in the ecliptic plane ; y-array 
normal to the sun , -y-array ojf-normal at an angle &i z an d Its 
A-frame at an angle &2 l a l A = 23.44°, orbit ascending node 
angle = 170 °, and inclination angle i = 285 ° 


i = o° 



OFF-NORMAL ANGLE 8 2z (DEG) 


Fig. 11. Annual roll-yaw momentum versus off-normal angle 
of -y-array; +y-array always normal to the sun; -170° 

5. Reaction Wheel Sizing and Configuration 
Optimization 

In Section 4, we observed that the momentum along the 
roll and yaw axes are coupled. Depending on the initial 
conditions, the secular momentum build-up tg a can exceed the 
capacity of the wheels either about the roll-axis, yaw-axis, or 
any other direction in the roll-yaw plane. Therefore, for the 
design purposes, the desired momentum capacity in the roll-yaw 
plane is the same in all directions. In contrast, the pitch 
momentum caused by the radiation torque is cyclic, with the 
amplitude equal to (gb/coo)- Hence, the pitch and roll-yaw 
momentum requirements leads to the well-known right circular 
cylinder momentum requirement for designing a wheel 
configuration, with cylinder axis along the pitch axis. On the 
other hand, the reaction wheels are also required to produce 
certain peak control torque about each of the three spacecraft 
axes. The torque requirements therefore form a rectangular 
parallelepiped. Both the momentum and the torque requirements 
can be met, in principle, using three wheels (not necessarily an 
orthogonal set); but for the sake of redundancy, four or more are 
often employed. In the following subsections, four-, six-, and 
three-wheel configurations are analyzed, two arrangements of the 
wheels and one-wheel failure are considered for each 
configuration. 

To visualize the wheel configuration most easily, first place 
the spin axes of all the wheels in the roll-yaw plane, perhaps 
radially symmetrically along the roll and yaw axes or otherwise, 
and then cant all spin axes, equally or unequally depending on 
the design, towards the pitch axis or its opposite. When the 
wheels are not along the spacecraft axes, a transformation 
matrix C bw is required to transform the wheel momentum 
vector H ww along the wheel axes to the total wheel momentum 
Hbw along the spacecraft axes: 

Hbw = Cbw Hww (^8) 
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When the number of wheels, n w , is more than 3, the matrix 
Cbw is rectangular, 3xn w , and its pseudo-inverse Cbw^ 

Cbw^ = Cbw^ (Cbw Cbw^) 1 (59) 

is required for the inverse transformation of (58): 

Hww = Cbw"*" Hbw (60) 

The above two transformation matrices are useful also for 
transforming the desired control torque Xc (3x1) about the 
spacecraft axes to the desired rate of change of the wheel angular 

momentum vector Hww (n w x 1): 


Hww = - £bw^ Ic (61) 

To determine the optimum cant angle(s) with the roll-yaw 
plane, it is logical to impose the requirement that, for a desired 
momentum vector capacity Hbw m spacecraft axes, the norm of 
the wheel momentum vector Hww be the least so as to 
minimize the cost and weight of the wheels. One suitable norm 
of the vector Hww is its Euclidean norm II Hww II defined by 


I Hww II = 


£h^ 


.i=i 


1/2 


(62) 


where H w j (i = 1 n w ) are the elements of the vector Hww- 

The minimization of II H ww II also results in minimum power 
consumption by the wheels for controlling the spacecraft, as 
shown below. 


Let co w i(0) be the initial speed of the wheel i, and 
<Dwi(t) the instantaneous speed while the wheel is acted upon by 
an electric motor, changing the wheel’s momentum H w i at the 

rate H w j(t). Then 


®wi(0 - ®wi(0) + Iw 1 j H w j (t) dt (63) 

where I w equals wheel’s moment of inertia about the spin-axis. 

For a constant H w j, the instantaneous power P w i consumed by 
the wheel i is given by, 

Pwi = I H w i I lo) w j(0) + I w _1 H w j 1 1 (64) 

where absolute values are taken to ensure that the power P w i is 

positive for both signs of H w j (t) and co w , (t).To determine the 
worst power consumption by the wheel assembly, we assume 

that C0wi(0) has the same sign as that of H w j. Furthermore, 
assuming that co w i(0) is the same for all wheels, the total 
power P w consumed by the n w wheels will be 


P w — I 0>w(0) I ^ I Hwi I + t Iw 1 HHww 11^ 

i = l 


(65) 


Thus, we see that, ignoring the linear term, minimization of 
power consumption leads to the minimization of the norm II 
Hww H. 

The vectors Hww and Hww are related, respectively, to the 
required momentum capacity H^w and the control torque 

capacity through the same pseudo- inverse matrix Q, w t. 

• 

Therefore, the minimization of both Hww and Hww yields the 
same optimum cam angle if the components of Hbw and Xc are 
proportional. This condition, however, is not always obeyed; 
for instance, in the presence of radiation torque, the desired roll 
and yaw momentum requirements are the same, but the desired 
roll and yaw torque requirements may be different because the 
roll and yaw errors might be controlled with controllers of 
different bandwidths, and the corresponding moments of inertia 
might be quite different. As a different example, the roll and 
yaw torques limits could be the same to facilitate momentum 
dumping with thrusters. The need for deliberation is thus 
evident. 


Four-Wheel Pyramid Configurations 
One possible arrangement of four wheels is shown in Fig. 12. 
The angle between two adjacent wheels is 90°, and they all are 
equally canted toward the -y axis by an angle q measured from 
the roll-yaw plane. When the cant angle q and the angle y in the 
roll-yaw plane are both zero, the momentum hi of the wheel 1 
is along the z-axis, h 2 along the x-axis, h 3 opposite to hj, and 
h 4 opposite to h 2 . The angle y is introduced so that the wheel 
torque can contribute, if desired, to all three axes and not just to 
roll and pitch or yaw and pitch. The corresponding rectangular 
matrix Qj W is 


£bw 


cqsy cqcy -cqsy -cqcy 
-sq -sq -sq -sq 

cqcy -cqsy -cqcy cqsy 


(66) 


whose pseudo-inverse, Eq. (59), is found to be 


sy/2cq — 1/4 sq cy/2cq 

r t cy/2cq -l/4sq -sy/2cq 

L0>w - -sy/2cq —1/4 sq -cy/2cq 

-cy/2cq -l/4sq sy/2cq 


(67) 


Let H x , H y , H z be the desired momentum capacity of the 
reaction wheels about the roll-, pitch-, and yaw-axis of the 
spacecraft. These three components disperse along the four 
wheel axes as follows, using Eqs. (60) and (67): 

"( H,sy + H z cy)/2cq - H y /4sq 
( H x cy - H z sy)/2cq - H y /4sq 

(~H x sy - H z cy)/2cq - Hy/4sq 

(-H x cy + H z sy)/2cq - H y /4sq 


Hww — 


( 68 ) 


Let, T x , Ty, T z be the desired maximum control torques about 
the x-, y-, and z-axis of the spacecraft. The maximum rate of 
change of the wheel angular momentum about wheel axes is 
then, according to Eqs. (61) and (67) 
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Fig. 12. Four-Wheel Configuration 


Hww 


(-T x sy - T z cy)/2cti + T y /4sr| 
(-T x cy + T 1 sy)/2cti + T y /4sTi 
( T x sy + T z cy)/2cn + T y /4sT| 
( T x cy - T z sy)/2cti + T y /4si) 


(69) 


If H x , Hy, H z are the desired momentum capacity of the 
reaction wheels about the roll, pitch, and yaw axis of the 
spacecraft, the 4x1 vector Hww is calculated analogously. Based 

on the definition (62), the norms II Hww H and II Hww H are then 
found to be: 


II H ww II 2 = (H x 2 + H z 2 ) / (20211) + Hy 2 / (4s2ti) (70) 


II Hww H 2 = (T x 2 + T z 2 ) / {2ch\) + T y 2 / (4s 2 !!) (71) 


which are independent of the angle y because the angle between 
two adjacent wheels is 90°. Comparing the two norms, it is 
clear that if (H x , Hy, and (T x , T y , are proportional, the 
minimization of one is the minimization of the other. The 
optimum cant angle q * is found to be 

tan 4 ii* = Hy 2 /2(H X 2 + H z 2 ) = T y 2 /2(T X 2 + T z 2 ) (72) 

which may be rewritten in a more revealing form: 

T x 2 / Ty 2 + T z 2 / Ty 2 = l/(2tan 4 il*) (73) 

which is the equation of a circle in the plane (T x /T y , 

T y ) and the radius of the circle equals 1/(^2 tan 2 q*). Eq. 
(73) states that as the torque requirements about x- and z-axis 
diminish, the radius of the circle shrinks and the optimum cant 
angle q* increases. This is exemplified in Fig. 13 where Eq. 
(73) is plotted for q* = 25°, 30° 60°. 


The norm II Hww 11 can be made dimension-free by dividing Eq. 
(71) with Ty 2 , This dimension- free right side of Eq. (71) is 
plotted in Fig. 14 as a function of q for a given (T x 2 + 
T z 2 )/T y 2 . The minimum value of the norm occurs at the 
optimum cant angle q*, and that minimum norm is found to be 

II Hww II 2 min / Ty 2 = [ ^2 (<T xy 2 + O zy 2 ) 1/2 + 1 ] 2 /4 (74) 

where the torque ratios a xy and are defined as 

G X y = T x / Ty G Z y = T Z / Ty (75) 

Table 2 furnishes optimum cant angle q* for several desired 
torque ratios; it also demonstrates that, for example, q* = 35.26 
as long as a xy 2 + a zy 2 = 2, regardless of the individual values 
of c x y and G zy . 

To determine the torque (or momentum) capacity of the 

wheels to produce the desired maximum torques (or 

• 

momentums) along the spacecraft axes, we again consider Hww, 
Eq. (69). Because T x , T y , T z are only three independent torque 

m 

requirements, the four elements of Hww are not all independent. 
Indeed, they are constrained by a relationship that is divulged 
from 


Hww — Cbw Cb w Hww- 


(76) 


(T x /Ty) 2 + (T Z /T y )2 = i/( 2 TAN4q) 


>-* 



Fig. 13. Dependence of minimum-power cant angle on desired 
torque ratios 
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Fig. 14. HtLwwil 2 1 T y 2 versus cam angle rjfor a 
four-wheel configuration 


Table 2. Examples of dependence of optimum cant angle on the 
desired torque ratios(<Jxy an ^ a iy ) 


HBfl 

■ 

1 



■zan 

0.5 

45° 

0.5 


1 

40.06° 

1 /V 2 

1 /V 2 

2 

35.26° 

1 

0 

<2 

1 

V2 

0 

4 

30.73° 

V2 

V 2 

8 

26.56° 

2 

2 


Recalling (66) and (67), Eq, (76) yields four identical equations 
staling 


H w i +H w 3 = H w 2 + H w 4 (77) 

The wheel momentums observe a similar relationship. 

Because of (77) and because of the nature of decomposition of 
the three desired torques T x , Ty, T z along the four wheel axes, 
Eq. (69), it may be intuitively clear that the magnitude of each 

element in (69) will not simultaneously reach H w max — the 
maximum achievable rate of change of wheel’s angular 
momentum. Therefore, recalling the definition (62) 

II Hww II 2 < 4H w , mx 2 (78) 

When the cam angle is optimum, the norm II Ji ww II is related 
to the specified torques T x , T y , T z according to Eq. (74). 
Therefore, in view of (78) 

[^(Tx 2 + T z 2)l/2 + T y ]/4<H w , mx (79) 


By way of illustration, if the torque requirements about the three 
axes arc all equal (T x = T y = T^, the inequality (79) yields 

Hw,mx > 3T x /4 (80) 

Considering that there are four wheels, each of capacity H w mx , 
for controlling three axes of the spacecraft each requiring the 
torque T x , the inequality (80) is perhaps a natural result, but it 

• 

does not reveal just how much H w mx must at last be, to size 
the wheel. For that, the Euclidean norm is not helpful and we 
must focus on Eq. (69) itself, as illustrated in the following two 
illustrations. 

Wheels Contributing to Roll and Pitch or Yaw and Pitch Only 

In Fig. 12, when y = 0 or 90°, each wheel contributes to 
either roll and pitch or yaw and pitch axes only. Regardless of 
y, the optimum cant angle for equal torque and momentum 
requirements satisfies, according to Eq. (72): 

T|* = 35.26°. Ian T]* = I/V 2 , sin r)* = iWT. cos T|* = V2/3 

(81) 

Substituting T x = T v = T z and st \* and CT|* from (81) in (69), 
we obtain, for y = 90° 

Hww = T x [—0.179 1.045 1.045 —0.179 ]^ (82a) 

which yields the desired maximum wheel torque capacity when 
not one wheel has failed 

Hw.mx - 1 .045 T x (82b) 

It is illuminating to compare the inequality (80) with (82b). 

Regarding the required momentum capacity of the wheels, 
we first obtain Hww, from Eq. (60), for y = 90°. Next, recall 
that from Eq. (41), the secular roll or yaw momentum at the 
end of one orbit is g a To (To = orbital period) and, from Eq. (42), 
the pitch momentum amplitude is gb/u> 0 - The least momentum 
capacity of a wheel for momentum dumping per orbit and for 
optimum cant angle (81) is then 

Hw.mx ^ ^3 Igalxo / 2 ^2 + V 3 lg b l / 4 cd 0 (83) 

Finally, in order to calculate the power consumption P w 

nw 

versus time t, Eq. (65), we require the quantities X 1 H wi I and il 

i=l 

Hww II 2 . From Eq. (82a), 
nw 

I I H wi I = 2.449 T x (84a) 

i=] 

which determines the intercept of the P w versus t curve at t = 0. 
The slope of this curve is proportional to the norm II H ww II 2 
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which, for equal torque requirement, is obtained from Eq. (74) 
by inserting a xy = 1 - a yz : 

II H ww II 2 = 2.25 T x 2 (84b) 

The indexes (84a) and (84b) should be kept as small as 
consistent with performance specifications. 

Wheels Contributing to All Three Axes 

In Fig. 12, when y = 45°, each of the four wheels 
contributes equally to the roll and yaw axes, as evident from Eq. 
(69). Reference 6 has examined this configuration to some 
depth. For y = 45°, for equal torque requirements (T x = T y = T z ) 
and optimum cant angle (81), Eq. (69) yields 

H ww = T x [ - V3/4 ^3/4 3 V3/4 V3/4 ] T (85a) 

The element with maximum absolute value yields the desired 
torque capacity of the wheel: 

Hw.mx * 3 V3 T mx /4 = 1.3 T x (85b) 

Comparing (85b) with (82), we conclude that, to produce a 
torque of magnitude T x about each of the three axes, the 
reaction wheels corresponding to y = 45° configuration must 
have 24.4% higher torque capacity than the wheels 
corresponding to y = 0° configuration. This is not surprising 
because for y = 45°, the wheel’s torque capacity is dispersed 
along all three axes, whereas in the case y = 0°, it is distributed 
along roll and pitch or yaw and pitch only. 

Following the derivation of the wheel momentum capacity, 
Eq. (83), the desired capacity for the configuration at hand is 


One-Wheel Failure 

For the four-wheel configuration shown in Fig. 12, we are 
usually interested in either y = 45° or 90° (y = 0° or 90° are 
effectively the same). And for these values, because all wheels 
are arranged symmetrically, failure of any wheel has the same 
consequences as the failure of any other. Therefore, to facilitate 
analysis, we arbitrarily assume the failure of wheel-3, and in 
that case the 3x4 transformation matrix Q, w , Eq. (66), 
condenses to a 3x3 matrix Cbw,3> formed by deleting the thud 
column of £bw-The inverse C.bw,3" 1 of Cbw,3 is found and 

then used to determine the torque vector H ww , following (61): 


Hww~ 


T x sy/crj - T z c y/cq 

T x (sy - cy)/2cr| + T y /2sq + T z (sy + cy)/2cq 
0 (wheel-3 failed) 

- T x (sy + cy)/2ct| + T y /2sq - T z (sy - cy)/2cr| - 


( 88 ) 


The wheel momentum vector Hww is determined likewise. For 
one-wheel failure case, the cant angle is not re-optimized 
because the cant angle of the wheels, once installed, is not 
changeable in the flight. 


Maximum Torque and Momentum Capacity When y~ 90° 

For y = 90°, and for the optimum cant angle (81), Eq. 
(88) yields 

Hww = T x [-1.225 2.091 0 0.866] T (89a) 


Hw.mx — ^3 ( lg a l IQ + Igb* / tOQ V 4 


(86) which in turn yields the required torque capacity of the wheel as 


where the optimum cant angle (81) has been used. Comparing 
(86) with (83), we find that H wmx now (y = 45°) is smaller 
than before (y = 0°), in contrast with the torque capacity 
conclusion drawn above. The reason of course is that the secular 
momentum toga is either about the roll-axis or yaw-axis, not 
both, whereas the torque capacity T x is desired about both roll 
and yaw axes. 

The two indexes of the power consumption are 
calculated with the aid of (85a) and (74): 


4 


I IH wi 1 = 2.598 T x 

(87a) 

i=l 


II Hww H 2 = 2.25 T x 2 

(87b) 


Comparing (87a) with (84a) and (87b) with (84b), we conclude 
that, for producing equal torque about the three spacecraft axes, 
and for the same initial wheel speed, the y = 45° wheel- 
configuration begins with a slightly higher power consumption 
and increases at the same rate as the y = 90° or 0° configuration. 


Hw.mx - 2.091 T x (89b) 

This is twice the required torque capacity in the no-failure case, 
Eq. (82b). 

Following the derivation of the momentum capacity Eq. 
(83) for the no-failure case, the momentum capacity for one- 
wheel failure case is: 


”V3 lg a l t 0 / ^2 

V3 lg a l T 0 / 2-il + \^3 lg b l /2o>o 
0 (wheel-3 failed) 

- V3 lg a l Xq/ 2^ 2 + V 3 lgbl /2(Oo - 


(90) 


Depending upon the relative magnitudes of g a and gb, either 
wheel- 1 or wheel-2 will yield the required momentum capacity 
(wheel-4 will yield the same capacity as the wheel-2). 

Regarding the two indexes of power consumption, Eq. (89a) 
furnishes 
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4 


X 1 H wi 1 = 4-182 T x 

(91a) 

i=l 


II Hww II 2 = 6.62 T x 2 

(91b) 


which may be compared with the no-failure results, Eqs. (87). 

Maximum Torque and Momentum Capacity When y- 45° 

For equal torque requirements about the roll, pitch, and yaw 
axes (T x = T y = T z ), and for the optimum cant angle r\* = 
35.26, the desired torque capacity of each wheel is 

Hw.mx - 1-732 T x (92) 

juxtaposed to the no-failure size (85b). By comparing the size 
(92) with the size (89b), the advantage of y = 45° configuration 
over Y = 0° or 90° configuration emerges: when one wheel fails, 
the y = 45° configuration can control the spacecraft with the 
wheels of smaller torque capacity than the y = 0° or 90° 
configuration can. The two indexes of power consumption are: 

4 

X I Hwi I = 3 V3T X = 5.196 T x 
i=l 

II Uww II 2 = 9 T x 2 (93b) 

Comparing (93) with (91), a disadvantage of the Y = 45° 
configuration is also unveiled: its power consumption is 
significantly greater than that of the y - 90° configuration. 
Finally, the desired momentum capacity is 

^w,mx ~ VI ( Iga 1 *0 + Igb 1 / “0) /2 (94) 

Compared with its no-failure counterpart, Eq. (86), the desired 
momentum capacity is now twice. 

Reference 6 may be reviewed for a different aspect regarding 
the selection of cant angle for the configuration at hand. 

Six-Wheel Pyramid Configurations 
Two-Cant-Angle Configuration 

One such configuration is shown in Fig. 15 where the 
wheels are arranged symmetrically (y = 60°), wheels 2 and 5 
controlling roll and pitch axes, and wheels 1, 3, 4, and 6 
controlling all three axes. Because of this fundamental difference 
between the two subsets of wheels, the cant angle r\2 of the 
former subset is allowed to be, in general, different from the 
cant angle qi of the latter subset. This freedom permits a 
greater economy in power consumption, if desired, and allows 
the reaction wheels to be of smaller torque and momentum 
capacity than the one-cant-angle configuration does. 

To determine the optimum cant angles and define 

Ci = cos Tii Si = sin T\ x (i = 1,2) (95) 

The transformation matrix (3x6) is 



~ c i/2 

c 2 

c i/2 

-ci/2 

-C2 -q/2 

Cbw~ 

-s l 

“ s 2 

1 

-s 1 

-si ~ s 1 


-V 3q/2 

0 

-V 3q/2 

-T 3cj/2 

0 V*3q/ 2_ 


(96) 


where, from the second and fifth column, it is apparent that the 
wheels 2 and 5 do not control the yaw axis, while the remaining 
four wheels control all axes. The pseudo-inverse matrix is 
determined using the definition (59): 



“ 

C c l/2> \ 


K ~*\ ^ 


r 1 ^ 

- 



c 2 


~*2 


0 


cj- 

< 

c ! /2 
-* i /2 

>Ac, 2 .2c 2 2 )^ 

-«1 

-•1 

i 2 + 2s 2 2 ) 

-1 

-1 

\T3cj 



_c 2 




0 



- 

^ c ]/2 J 







(97) 

and the Euclidean norm of the vector H ww is 


H Hww H 2 - T x 2 /(ci 2 + 2c2 2 ) + 

T y 2 /2(2si 2 + S2 2 ) + T z 2 /3ci 2 (98) 

which is minimized by the optimum angles r|i* and T)2* 
defined by 

5^1*= (Tx + Ty - T 2 ) / (T x + T y + T z ), 

c^Tjl* = 2T z /(T x + T y + T z ) (99a) 

5^2*= {2(T Z - T x ) + Ty} / (T x + Ty + T z ), 

ch\ 2 * = (3T X - / (T x + T y + Tz) (99b) 

Fig. 16 portrays the optimum angles r\\* and rj2* for specified 
torque ratios a X2 and a yz . Substituting the optimum 
trigonometric functions in Eq.(98), the minimum value of the 

norm II Hww H 2 is found to be 

II Hww II 2 = (T x + T y + Tz) 2 /6 (100) 

The three torque components T x , T y , T z are independent, and 

once specified, they are produced by the six Hwi (i=l, 6), 

given by Eq. (61). Clearly, these six quantities are constrained 
by three relations which are obtained from the expanded version 
of Eq. (76). Due to these constraints, 

II Hww II 2 < 6 H w mx 2 (101) 

analogous to the inequality (78) for the 4-wheel configurations. 
Combining (101) with (100), we obtain 

(T x + T y + T^ / 6 < H w mx (102) 
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Hw,mx > T x /2 


(103) 



Fig. 15. Six-wheel hexagonal configuration 



m*«o° 

(Oyz “ — C?XZ+1) 


Fig. 16. Dependence of optimum cant angles Tjj* and Tf2 * on 
the torque ratios a u and a yz 

which states that the sum of the maximum torques that can be 
produced about the three spacecraft axes must be less than the 
total torque capacity of the wheels — the cant angles being the 
underlying reason. For equal torque requirements (T x = T y = 
Tz), the inequality (102) reduces to 


which may be compared with (80). It is intuitively clear that, 
instead of arranging six wheels as shown in Fig. 15, if they 
were arranged two wheels per axis, then for equal torque 
requirement about the three axes each wheel’s torque capacity 
must satisfy 

H w> mx - T x /2 (104) 

instead of (103). 

When the maximum required torque components T x , T y , T z 
are all equal, Eqs. (99) yield Eq. (81) 

STij* = l />/3 = sti 2 *: ct) ] * = V2/V3 =cn 2 *; 

Tl t * = 35.26° = n 2 * (105) 

The two cant angles, therefore, coalesce and indeed they become 
the same as that for the 4-wheel configurations. The desired 
maximum torque capacity of the wheel is then found to be 

H w ,mx ^ 0.846 T x (106) 

which is smaller than the torque capacity (82b) or (85b) for 4- 
wheel configurations for the same torque requirements about the 
spacecraft axes. 

Regarding the power consumption, the Euclidean norm 
(100) yields 

IIHwwII 2 = 1-5 T x 2 (107) 

Comparing (107) with (84b) and (87b), we observe that for the 
same T x about all three axes, the power consumption of the 6- 
wheel configuration increases at a smaller rate than that of the 
4-wheel configuration. Finally, 

6 

X I H w i I = 2.509 T x (108) 

i=l 

which is within the two values (84a) and (87a) for the two 4- 
wheel configurations. 

One-Cant-Angle Configuration 

A hexagonal wheel assembly with two different cant angles 
might be difficult to install in a spacecraft bus; so we now 
optimize a hexagonal configuration with one cant angle. The 
pseudo-inverse matrix, Eq. (97), simplifies and Eq. (61) yields 
the Euclidean norm 

II H ww II 2 = (T x 2 + T z 2 ) / Sc 2 !! + T y 2 /6s 2 ti (109) 

where the subscript 1 of T|i is dropped because now there is 
only one cant angle. Minimization of this norm leads to, 
surprisingly, the condition (72) for the 4-wheel configurations. 
The minimum value of the norm (109) is 

11 Hww II 2 = [Ty + V2 (T x 2 + T z 2 ) ] 2 / 6 (110) 
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Fig . 17. A Three-Wheel Pyramid Configuration 

which is two-thirds of the value in (74) for the 4-wheel 
configuration. Also, because of the inequality (101), we arrive 
at 


Hw.mx > [V(T X 2 + T z 2 )/18 + T y / 6] (1 1 1) 

which is different form (102) but, for equal torque requirement, 
reduces to (103). 

The one-wheel failure analysis, considered above for the 
four-wheel configurations, becomes unwieldy because of the 
5x3 size of the reduced matrix £b w . Therefore, pertinent results 
such as wheel torque capacity and power consumption are 
obtained with the aid of a computer and summarized later in 
Table 3. 

Three-Wheel Pyramid Configuration 

When wheel redundancy is not warranted, when for reasons 
of cost and weight the number of wheels must be bare 
minimum, and when the torque requirements about the three 
axes are not necessarily equal, the three- wheel pyramid 
configuration shown in Fig. 17 might be an ideal choice. For 
this arrangement, the three wheel angular momentums can be 
expressed in spacecraft axes as follows: 


-H* - 


-CT| 1/2 C T| 1/2 CT] 


"H wl - 

Hy 


-ST| -ST| -ST| 


H*2 

-H* - 


0 -V~3/2 CTj \^ 3/2 CTJ 


-h w3 - 


wherein the 3x3 transformation matrix is £bw- Because now 
there is no redundancy, the pseudo-inverse matrix £Lb w ^ 
becomes the regular inverse matrix £b w -1 - After determining 


, • 

C.bw'' 1 and substituting that in Eq. (61), the vector Hww in 

terms of the required torque components T x , T y , T z turns out to 
be 


Hu 


“2T x /3cr| + T y /3sTi 
-T x /3cr| + T y /3sq + T z /3cq 
_-T x /3cq + T y /3sq - T z /3cq . 


(H3) 


For minimization of power consumption, we arrive at the 
followin Euclidean norm. 


II Hww II 2 = 2 (T x 2 + T z 2 ) / 3c 2 !) + Ty 2 ^ (1 14) 

which is four-thirds of the norm (71) for 4- wheel configurations 
and six-thirds (twice) of the norm (109) for 6-wheel 
configuration, implying that if they all begin from zero wheel 
speed, the 3-wheel configuration will consume greater power in 
the stated ratio. For example, for equal torque requirements (T* = 
T y = T z ), while the minimum value of the 4-wheel norm is 
2.25, that of the 3-wheel norm is 3.0, which is, incidentally, 
the same as that for the three orthogonal wheels one per axis. 
Moreover, if the cant angle is not set to be the optimum (rj * 
q*), the three-wheel pyramid configuration will use more power 
than the one-per-axis configuration. Next, the norm (114) yields 
the same optimum angle as one for the 4- and 6-wheel 
configurations, Eq. (72). For this optimum angle, the following 
minimum value of the norm emerges: 

N Hww min = [ T y + V 2(T X ^ + T 2 ^) ] /3 ( 115 ) 

which is twice the value (1 10) for the 6- wheel configuration and 
four-thirds of the value (71) for the 4-wheel configuration. 

The required torque capacity of the wheels for equal torque 
requirements about roll, pitch, and yaw axes is found to be 

H w ,mx > 1-39 T x (116) 


As a check, note that the norm of Hww for three-wheel 
configuration is indeed 


II H ww II 2 = 3T X 2 (117) 

equal to that for a one-wheel-per-axis configuration; but the 
power consumption of each wheel would be quite different from 
that for the one-wheel-per-axis configuration. 

Overall Comparison of Six Configurations 

When the torque requirements about the roll, pitch, and yaw 
axes are not the same, the wheels of different torque capacities 
along different axes might be selected; but from the standpoint 
of reliability and cost, that is usually not preferred. Perhaps 
a more attractive choice is a six -wheel configuration with 
identical wheels, the cant angle selected according to the desired 
torque ratios. For equal torque requirements, the optimum cant 
angle is q* = 35.26°, and the associated wheel torque capacity 
for the required torque T mx must be at least 0.846 T mx [Eq. 
(106)] — greater than 0.5 T mx for the two-wheel-per-axis 
arrangement The two power consumption indexes in the case of 
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no wheel failure shown in the second row of Table 3 restate 
Eqs. (107) and (108). When the wheel in the roll-pitch plane 
fails, the torque capacity of the remaining five wheels must be 
boosted to at least 1.31 1 T mx to produce the required torque T mx 
about the spacecraft axes. This result is obtained by failing the 

wheels 1,2 6, one at a time, and then determining the 

absolute maximum value of the wheel torque in each case for 
generating T m * torque about each of the three spacecraft axes. 
The maximum Euclidean norm and the associated absolute sum 
6 

£lH w ji are also shown in Table 3. Comparing the 2-wheel-per- 
i=l 

axis and 6-wheel hexagon configurations, we find that for equal 
torque requirements , the latter (hexagon) configuration requires 
wheels of larger torque capacity and it consumes more power — 
and therefore not as favored as— the former configuration. 
However, when the roll, pitch, and yaw torque requirements are 
not the same, the conclusion will possibly swing in favor of 
the hexagon configuration. 

Although six-wheel configurations provide substantial 
reliability and three-wheel redundancy, they could be expensive, 
so four-wheel configurations may be desirable instead, which 
provide a one-wheel redundancy. Two such configurations — one 
with pyramid base parallel to roll-yaw axes and the other with 
the base at 45° — are discussed above. For the purpose of 
comparison , call these configurations parallel- and 45 °- 
configuration, respectively. Under no-failure case, the 45°- 
configuration requires wheels of larger torque capacity than the 
parallel -configuration, but in the event of a one wheel failure, 
the situation reverses. On the other hand, from the power 
consumption viewpoint, under no-failure case, the 45°- 
configuration uses only slightly more power than the parallel- 
configuration, but the failure of a wheel aggravates this 
difference. Because the final design is usually based on one- 
wheel failure performance, we infer that if power is relatively 
abundant and the wheel torque capacity is at a premium, the 45°- 


configuration should be selected. On the other hand, if power is 
expensive and the cost of the wheels depends only weakly on its 
torque capacity, the parallel configuration will then be a more 
prudent choice. 

When wheel redundancy is not warranted, only three 
wheels — necessary and sufficient for spacecraft control — can be 
employed. If the torque and momentum requirements about the 
three axes are identical, the control engineer may opt for one- 
wheel-per-axis configuration. But in the case of dissimilar 
requirements, 3-wheel pyramid, with the cant angle suitable to 
the desired torque and momentum ratios, might be preferred. 
Table 3 compares these two 3-wheel configurations for equal 
torque requirements, and shows that the pyramid configuration 
requires wheels of 39% bigger torque capacity, although its 
power consumption may be slightly less than or equal to that of 
the one-wheel-per-axis configuration. 

Fig. 18 sums up the comparison between the power 
consumption of the six configurations considered in Table 3 for 
equal torque requirements. In particular, the Euclidean norm of 

the vector H ww versus the cant angle for each configuration for 
the no-wheel-failure case are shown in the figure. As noted 
before, the cant angle q* (q* = 35.26°) for minimum power 
consumption is the same for 3-, 4-, or 6-wheel pyramid 
configurations. 

Concluding Remarks 

Among a variety of disturbance torque that act on a 
space-craft, only solar radiation is considered in the preceding. 
For clarity, the torque expressions are further specialized by 
assuming that the vehicle mass center always remains in the 
pitch-yaw plane. Although this was true for the spacecraft that 
led to this study, the roll component of the vector from 
instantaneous vehicle mass center to the geometric center of the 
array or bus may not be zero for other spacecraft. Also, while 
solar torque vanes at orbit frequency, aerodynamic torque, for 


Table 3. Comparison of Six Configurations for Equal Torque Requirements about Roll, Pitch , and Yaw Axes, 

Based on Minimum Power Consumption 

REQUIRED CONTROL TORQUE: T cx = T cy = T c2 - T m * 


CONFIGURATION 

OPTIMUM 

ANGLE 

NO-FAILURE 

WORST 1 WHEEL FAILURE 

TORQUE 

POWER 


Required 

Torque 

Capacity 

Hw-mx^mx 

Total Power 
Consumption Rate 
r>w 

hi 

Total Power Intercept 
Due to Nonzero Initial 
Wheel Speed 
r>w-l 

Ifyl^mx 

tel 

Faied 

Wteel 

I 

WwH 

Producing 

Max 

Torque 

Requrefl Torque 
Capacity 

^Vm/Tmx 

Fated 

Wheel 

* 

Total Power 
Consumption 
Rate 

rw 

X ^ 2/t ™ 2 

tel 

Total Power Intercept 
Due to Nonzero Initial 
Wheel Speed 
rv-i 

l^il/Tm* 

tel 

2 Wheel/ Axis 
("W*J 


0.5 

15 

3.0 

Any One 


1.0 

Any 

One 

2.0 

3 

6- Wheel Hexagon 
(rv-6) 

35.26 

0.846 

15 

2.509 

#5 

#4 

1.311 

# 4 

2.933 

3073 

4- Wheel Pyramid. 
Base Edges II to 
Axes (n*=4) 

35.26 

1.045 

2.25 

2.449 

#1 

#2 
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Fig. 18. Reaction wheel configurations trade-offs for equal roll, 
pitch, and yaw torque requirements 

instance, may vary at twice the orbit frequency and a yaw bias 
torque might arise. For these different circumstances, the torque 
and momentum expressions have to be derived afresh to size the 
wheels. Regarding the wheel configurations, besides the two 
four-wheel configurations considered in the paper, there are two 
more: 1) NASA's standard four-wheel arrangement of one wheel 
along each body-axis and the fourth wheel inclined equally to all 
three axes; and 2) all four wheels canted equally to the pitch axis 
and each controlling roll and yaw as well, but more inclined to 
roll-axis than to yaw -axis or the converse depending on the roll 
and yaw unequal torque requirements. In the first arrangement, 
the cant angle of the fourth wheel is already determined, only 
the wheels' torque and momentum capacity need to be sized for 
one-wheel failure scenario. In contrast, in the second 


arrangement, two angles must be optimized to minimize power 
consumption for given torque and momentum requirements: the 
cant angle q with the roll-yaw plane and the angle y with the 
roll axis for all four wheels. For these optimum angles, the 
torque and momentum capacity of the wheels will be sized 
according to one- wheel failure condition, as shown in the paper. 
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